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Abstract 

We present the Mathematica application .DoFwrl^which allows to derive Dyson-Schwinger equations 
and renormalization group flow equations for n-point functions in a simple manner. DoFun offers 
several tools which considerably simplify the derivation of these equations from a given physical 
action. We discuss the application of DoFun by means of two different types of quantum field 
theories, namely a bosonic 0{N) theory and the Gross-Neveu model. 

Keywords: Dyson-Schwinger equations, functional renormalization group equations, correlation 
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PROGRAM SUMMARY 

Program Title: DoFun 
Version number: 2.0.0 

Licensing provisions: CPC non-profit use license 
Programming language: Mathematica 7 and higher 

Operating system: all on which Mathematica is available (Windows, Unix, MacOS) 
PACS: 11.10.-z,03.70.-fk,11.15.Tk 

Nature of problem: Derivation of functional renormalization group equations and Dyson-Schwinger equa- 
tions from the action of a given theory. 

Unusual features: The results can be plotted as Feynman diagrams in Mathematica. The output is compat- 
ible with the syntax of many other programs and is therefore suitable for further (algebraic) computations. 



1. Introduction 

The derivation of Dyson-Schwinger equations (DSEs) around 1950 [HE] and renormalization 
group equations (RGEs) in the early 1970s gHS] has equipped us with powerful tools for an analysis 
of the dynamics of quantum field theories. Both approaches have been further developed in the 
past 30 years. In fact, many formulations of these two methods now rely on a formulation in terms 
of so-called generating functionals for Green functions [101414) . 

These days functional approaches, such as DSEs, RGEs or the n-PI formalism, see e. g. [151 - 
[T5] , are well-established for studies of quantum field theories. Apart from functional approaches, 
Monte-Carlo simulations based on a discretized action have been extensively used to study non- 
perturbative phenomena. In fact, so-called QCD lattice simulations are currently the most powerful 
tool available for a study of full QCD. However, the implementation of fermions in such simulations 
continues to be a non-trivial task. Functional approaches are also non-perturbativc but do not 
have problems arising from a discretized action or from the implementation of fermionic degrees 
of freedom. However, the application of DSEs and RGEs eventually requires in most cases a 
truncation of the full system of equations of a given theory. From this point of view, it is clear that 
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Monte-Carlo simulations and functional approaches are complementary approaches for studies of 
non-perturbativc phenomena in quantum field theories. 

DSEs and non-perturbative RGEs have been successfully employed to gain a better under- 
standing of a large and diverse variety of quantum field theories. For instance, detailed studies 
of condensed-matter systems, see e. g. [T1H2I], critical phenomena, see e. g. [I^HIS], few- and 
many-body physics, see e. g. p7H52] . gravity, see e. g. [55H57] . QCD, see e. g. |55H75] . standard 
model physics, see e. g. |75J|77], and supersymmetry, see e. g. [75H5T] . are available these days. 
Furthermore, DSEs are also used as an alternative to the traditional Feynman graph approach in 
perturbation theory, see for example |82H84j . For reviews on and introductions to the application 
of DSEs and non-perturbative RGEs we refer to e. g. gH [Ml iSl [86] and [l6 t [87H96 ] . respectively. 

The applicability of DSEs and RGEs to very different theories is indeed an attractive feature of 
these approaches. However, the intricacy of the derivation of DSEs and RGEs scales non-linearly 
with the complexity of the theory. Therefore pushing a computation to a higher level of accuracy 
often requires a big effort as the number of terms increases considerably. Consequently, tools 
facilitating the derivation of such equations are helpful for future studies with DSEs and non- 
perturbative RGEs and our Mathematica 97J program _Dofw7|^does exactly that: It allows for an 
automatic derivation of DSEs and non-perturbative RGEs from a given action. Of course, finding 
a suitable ansatz for the effective action for a given problem is left to the user and remains to be 
the most difficult step from a physical point of view. 

The program DoFun is a further development of the Mathematica DoDSE package]^ [98] which 
was restricted to a derivation of DSEs in symbolic form. There are two new main features: First, 
the derivation of RGEs is now included and, second, the symbolic results can be transformed into 
the corresponding algebraic expressions. For the latter task two additional packages were added 
to DoDSE: DoAE and DoFiQ To account for the inclusion of RGEs the former package DoDSE 
was renamed DoDSERGE. These three packages form the content of the program DoFun. 

In the process of the development of DoFun we always had in mind that there indeed exists a 
variety of programs to deal with single steps of the derivation of functional equations like performing 
Dirac traces or simplifications of the color structure. We did not intend to force the user to learn 
an additional new program syntax for these kinds of tasks but to open up the possibility of using 
our program in combination with the corresponding available programs. Therefore, our goal was to 
stay as general as possible and to allow for the combination of DoFun with many well-established 
programs, such as TRACER or FeynCalc |100j. The program DoFun only performs the most 
basic simplifications and the output can then be handled with other programs. For example, the 
user may want to stick to his favorite program to deal with the color algebra or the user may want 
to use his own programs for certain operations. The latter is often very efficient, since the code 
can then be specifically tailored for the problem at hand. Following this general approach allows 
a high flexibility for the user and the treatment of a huge variety of theories. 

We are not aware of the existence of other programs for the derivation of DSEs whereas other 
packages for the derivation of functional RGEs indeed exist, see Refs. jl01H103 . However, as men- 
tioned above, DoFun is based on the application DoDSE and exploits similarities in the derivation 
of DSEs and RGEs. It thus can be viewed as complementary to other existing programs. 

The aim of DoFun is not to replace the manual derivation process completely, but to provide 
additional help in cases where the manual derivation becomes too cumbersome. While 0{N) models 
may still be more easily accessible with pen and paper, studies of more complicated theories may 
benefit in various ways from DoFun: 

• In general, a derivation of functional equations by hand becomes tedious when several differ- 
ent fields or large numbers of operators are taken into account. 

• With a package for an automatized derivation of functional equations, equations at the next 
higher level of the truncation are easily accessible. For example, this is particularly useful to 
explore the role of higher-order operators. 

• Theories with complicated tensor structures benefit from the connection to a computer alge- 
bra system as tensors can be directly computed or simplified. 



■^Short for "Derivation of Functional equations". 

^ Short for "Derivation of Dyson-Schwinger Equations". 

*The names are abbreviations of "Derivation of Algebraic Expressions" and "Derivation of Feynman Rules", 
respectively. 
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Figure 1: Workflow of DoFun. 



• Finally, the graphical output of DoFun is often helpful for an illustration of the basic structure 
of the equations. 

DoFun is an extension of DoDSE. Both have already proven to be useful for a number of 
studies. In fact, DoDSE was first used for the derivation of the DSEs of the maximally Abelian 
gauge [Ml I104| . Actually, the complexity of these equations was the reason for the development 
of DoDSE in the first place. While a manual derivation of the equations would still have been 
possible but very ineffective, the computation of certain quantities would have been out of reach 
without the aid of a computer algebra system |105j . The advantages of an automatized derivation 
of equations has also been appreciated in the analysis of the Gribov-Zwanziger action, where one 
has to handle large expressions [SSI 11041 1106j . Studies of scalar fields coupled to Yang-Mills theory 
benefited from DoDSE as well ' 107H110] . Recently, DoFun has been used for an investigation of 
bound states appearing in the BRST quartets of QCD [1111 1112j . 

Working with DoFun involves several steps: The first one is the derivation of the equations in 
symbolic form using the functions doDSE or doRGE. As input they require the action in symbolic 
form. In order to transform the results into algebraic expressions with DoAE, Feynman rules have 
to be defined. They can be either derived by hand from the physical action or with the aid of other 
available packages, such as DoFR. The algebraic expressions represent the final output of DoFun. 
Further manipulations are up to the user. The generic workflow with DoFun is summarized in 

To make DoFun easily accessible for the user we implemented documentation directly into 
the Mathematica help system. For quick reminders of the syntax of a function one can use the 
command ?, e. g. 

?doRGE 

which also provides simple examples for many functions. For more detailed information we also 
included a section into the Documentation Center of Mathematica. There, tutorials and guides 
introducing the basic features of DoFun can be found. It can be accessed via Help — > Documentation 
Center — > Add- 0ns and Packages — > DoFun -> Documentation. All functions of DoFun have a 
help entry which can be opened by moving the cursor into or behind the function name and pressing 
the key Fl. 

To install DoFun the directory DoFun should be copied into the Mathematica subdirectory 
Applicationsj^ Within Mathematica the program DoFun can then be loaded with <<DoFuri or 
Needs ["DoFuri"] and the help system is accessible. 

The paper is organized as follows: In Sect. [2| we briefly summarize the derivation of DSEs and 
non-perturbative (functional) RGEs. This part can be skipped by readers familiar with the topic. 
In Sect. [3|the basic constructs and quantities of the program DoFun are introduced. In Sect. [4| 
we explain the main functions of DoFun. Their usage is then demonstrated in Sect. [5|by deriving 
the flow equations of two substantially different theories, namely an 0{N) model in d space-time 
dimensions and the Gross-Neveu model in = 3 space-time dimensions. Comments concerning 
the implementation of fermions are given in [Appendix A[ In [Appendix B[ we give a short version 
summary of DoFun. 



^Under a Unix system tliis is typically 7 .Mathematica/Applications. 
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2. EXinctional RG equations and Dyson- Schwinger equations 



2.1. Actions and correlation Junctions 

The program DoFun deals with the effective action r[$] as the underlying quantity for the 
derivation of correlation functions. The field $ denotes a collective field vector, e. g. $ = 
{(fi^(x),...,Lp^{x)} for an 0{N) model or $ = c°(a;), c°(x), g^(a;)} for QCD. Here 

and in the following we shall assume that the index of the field $ refers to the field type, the 
momentum and to all indices associated with internal symmetry groups (except if these indices are 
stated explicitly). The effective action is defined as follows: 

r[$] := sup(-Ty[J] + J,) , (1) 
J 

where the J^'s denote sources for the fields $i. Note that we sum/integrate over all possible values 
of index variables that appear twice in a single term. On the other hand, the generating functional 
W[J] is related to the bare action S'[(/)] via the path integral: 

Z[J] = J D[(j)]e-'^+'f''''' =: (2) 

Here, the ^^'s denote the quantum fields. The so-called average field $ is given by 

SW 



Z[J]-^ j D[<i>]c^,e-''+^^J^ . (3) 



The physical expectation value of the fields 4>i is then obtained by setting the external source J to 
zero, $phys := {(l>i)j=o- 

The effective action can be expanded about, e. g., the physical ground state <&phys as follows: 

oo ^ 

rt'J'] = E Jf^^ E r'"-""(*n - ^.upi^ys) ■ • . (<f^,. - $.„,phys) , (4) 
n— 

where jV*i - '" is the corresponding symmetry factor and the (physical) n-point functions are given 
by the expansion coefficient^ 



(5a) 

^ — ^phys 



(5b) 

* = *Dhvs 



Note that the physical ground-state may even be space-time dependent as recently discussed in 
the context of Gross-Neveu and Nambu-Jona-Lasinio models, see e. g. Refs. |113H116] . 

From the generating functional W[J] of connected correlation functions we obtain the propa- 
gator, i. e. the inverse of the two-point function: 



{^iy' > (6) 



where the index J denotes the dependence of the two-point function on the source jj^ From this 
it follows that the derivative of the propagator with respect to a field yields two propagators and 



®For convenience, we have inserted a minus sign appearing in our definition of the vertices T'l "'". With this 
definition, the sign of any diagram can be obtained straightforwardly as it does not depend on the number of vertices 
in a diagram. 

""Note that (5r/<5$i = J;. 
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a three-point functiorj^ Overall, we only require the following set of differentiation rules: 
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^3 = 5^J, (7b) 
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These relations form our basis for an algorithm for the derivation of DSEs and functional RGEs 
which underlies DoFun. 

2.2. Derivation of Junctional renormalization group equations 

Let us now briefly discuss the derivation of functional RG flow equations. Here, we shall follow 
the standard derivation of the flow equation for the so-called effective average action given in 

Ref. ng. 

First, we define the so-called effective average action which depends on a momentum-shell 
parameter k and interpolates between the bare (classical) action and the full quantum action F 
for fc — >■ 0. To this end, we introduce a so-called cutoff function R}^ into the path integral ^ as 
follows: 

^ + ASk [0] = + . (8) 

The so introduced (momentum) scale k allows us to integrate out quantum fluctuations in a 
controlled way. In momentum space the cutoff action ASk[(j>] can be written as follows 

1 , If d'^q d'^q' 



2 



^Rfiq)cba{q)M-q)- 



Note that the regulator function Rk is matrix-valued in field-space and has to obey the following 
constraints: 

(1) limg2/j,2^o^fc > 0, 

(2) limfc2/^2^o^fc =0, 

(3) limfc^A^oo i?fc = oo . 

The first constraint implements an IR regularization for the path integral. The second constraint 
ensures that the regulator vanishes for fc -> 0. Thus, we recover the standard generating functional 
Z defined in Eq. ^ ioi k 0. 

The inclusion of the cutoff action ASk renders all generating functionals fc-dependent. In partic- 
ular, the effective action now depends on the momentum scale fc ~ q. The scale-dependent effective 
action F^, the so-called effective average action, is defined via a modified Legendre transformation: 

Ffc[$] = -Wk[J] + J^<^>^ - l^^R]^^J (10) 



with 



(11) 



Note that (5(MAf-i)/5<I> = 0. 
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By construction, the effective average action Ffc includes all quantum fluctuations associated with 
momenta p ^ k. The modified Legendre transformation together with the third constraint for the 
regulator function ensures that reduces to the classical action 5 for fc — > A. On the other hand, 
we obtain the full quantum effective action F in the limit fc — > 0. 

From the effective average action F^ we may derive scale-dependent correlation functions " 
by differentiating F^ with respect to the fields]^ Thus, the differentiation rules given in Eq. Q 
also apply to the effective average action Ffc. 



Taking the derivative of Eq. ( 10 1 with respect to the scale k we obtain the RG flow equation 



for the effective action, the so-called Wetterich equation \T2\ : 

dkTkm =\ [{Vkm+Rk)-^Y ^^^^ 

= isTr(Ffc[$]+i?,)"'4i?fc. (12) 

Clearly, the function i?^, implements an IR regularization for the momentum integrals and specifies 
the details of the Wilsonian momentum-shell integrations. Note that the super trace STr includes 
a minus sign for Grassmann- valued fields. 

The package DoDSERGE for the derivation of the flow equations of ?i-point functions is based 



on a simple reformulation of Eq. ( 12 ) 



atr,[$]=^STr4in(r,[$] + i?fe), (13) 

where t — ln(fc/A) with A being a UV cutoff scale. The derivative dt only acts on the regulator 
Rk ■ This formulation of the flow equation represents the starting point for the derivation of RG 
flow equations with DoFun. The derivative dt is then taken at the end of the generation of the flow 
equations. This last step increases the number of diagrams since it generates diagrams with the 
same topology but with the regulator insertions att ached to different internal lines. Note that it is 



possible to suppress the derivative dt, see Sect. 4.2 This is in fact useful if we are only interested 
in checking the structure of an equation. 

The flow equations for, e. g., mass terms or so-called wave-function renormalizations are 
extracted from the flow equations of the n-point correlation functions. To obtain them, we take 



functional derivatives of Eq. ( 13 1 with respect to the flelds. The first derivative yields 



^rdi^^m =;^STr Ut ,J]l'^ D'A . (14) 

where the indices i and j are contracted by the trace operation. From this, higher derivatives are 
then obtained by using the differentiation rules Q. By definition, the derivative dt only acts on 
the propagators and inserts the derivative of a regulator into the loops: 

dtD'^ = -D'^idtR^^D""^. (15) 

By taking derivatives with respect to the fields, the flow equations for the n-point functions 
can be derived in an exact form yielding an inflnite tower of coupled flow equations. Since it is 
impossible to study the flow of an effective action containing all operators allowed by the symmetries 
of a given theory, we have to restrict ourselves to that subspace of operators which we expect to 
be (most) relevant for the physical problem under consideration. This is the most difficult step 
since it requires a lot of physical insight into the problem in order to choose the correct subspace 
of operators. 

We now turn to a discussion of the command doRGE which performs the derivation of functional 
RGEs. It follows these steps: 



(1) Take the first derivative of dtTk, see Eq. (14) 



(2) Take higher derivatives using the command derivRGE. Note that this command takes into 
account the proper ordering of fermion fields. 

(3) The sources are set to their physical values with setSourcesZeroRGE. 



Prom here on we suppress the index J as it should be clear from the context if an expression depends on J 7^ 0. 
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(4) Take the derivative with respect to t ^ ln(fc/A), see Eq. (15) 



(5) Identical Feynman diagrams are summed up with identif yGraphsRGE. 

Let us now give an exphcit example of the steps performed by the command doRGE. To this 
end, we derive the flow equation for the ghost two-point function in Yang-Mills theory in the 
Landau gauge. There we have two types of fields, namely gluons and ghosts denoted by A and c. 
The collective field vector $ is given by $ = {A°,c°,c"}. Note that the ghost fields c and c are 
Grassmann- valued, i. e. anti-commuting, fields. The extra rules required for their treatment are 
described in [Appendix A To obtain the ghost two-point function we first take the derivative of 



Eq. 13 with respect to a ghost and an anti-ghost field. This yields (step 1) 



5Ca Scb 2 6Ca ^ ^ /J Scb 



STr^t ,, A.,—. (16) 



The indices a and h are collective indices including the momenta and the color indices of the 
fields. On the other hand, the indices «, n and j of the collective field $ also include a field index 
which refers to the type of the field (gauge field, ghost, anti-ghost). The open indices i and j are 
contracted by the trace operation but we perform the trace only at the end when we know if the 
corresponding fields are Grassmannian or not. Performing the second derivative we obtain (step 

2): 

5Ca 6cb 2 \\ 6CaS<^,6^n 5cb j 5Ca5^i5<^n \ Scb 

Before taking the (super-) trace we set the external sources to zero (step 3): 

1 , x4p f4p r4p 
d pa*>_ — -STr9 ■! " ^ fc jjnj _j " fc j^nj _| " ^ fc jjnj 

^''^'^ 2 \5c.n5Ai5A„5cb Sc„5ci5cnSch '^'^ 5c„6cj5c„5ch ^'^ 



^^^k j^nm '^'^^fc T^l] ^^^k „nm ^^^k „ij \ 



SCaSAiSCn SCmSAiScb SCaSciSAn ^^SA^Sci6cb 



^ k.cAAc^ AA k.cccc-'-^cc ^ k,cccc^cc 

T^k,cAc i-^nm-r^mlb r-dj ^ T^k,ccA p.nm-pmlb r)0 ' 

ain cc ^ k,cAc-^AA ain AA k.Acc cc \ 



isTrO.{ 

}• (18) 



2 I k,cAAc^AA ^ k,cccc^cc k,cccc^cc 



-pk.cAc j^nni-pralh j^lj i -pk^ccA j^nm-pmlb T)3^ 

^ ain ^cc ^ k,cAc^ AA ' ^ ain ^ AA'- k,Acc^cc. i 

In the last step we have brought the ghosts into canonical order. By setting i — j we perform the 
trace. However, we have to take into account a minus sign if i and j correspond to indices of ghost 
fields: 

dtTtc = ( - ^'^IacDTa + ^TI^Ld^: - 2ri^J,,D-^Ti%^D^J,^) . (19) 



By summing up the terms in Eq. ( 18 ) we anticipated step 5. Finally we have to take the derivative 



of the propagators with respect to the RG time t (step 4): 

dj^cc 1 pamfc r^nt /o rylni \ Tinii painb T^nl ( c\ rtlrn \ r^mi 

k,ab — 2 k,cAAc^ AA\"tnf. j^j^jU ji^js^ — i k,cccc^cc\(^t^k,cc)^cc 

.-pian T^nni-plmb j-^lj / jjjo \ 7-)°* , -pian riTij / c, jyjo \ -r^om-plmb 7-)H ( 0(W 

+ ^ k,Acc^cc ^k,Acc^AA\"t^k,AA>^AA + ^k,Acc^cc\Ct^k,cc)^cc ^ k.Acc^ AA- l^^lj 

This is the flow equation for the ghost two-point function. It is depicted in Fig. [2j We observe 
that the pure ghost loop differs by a minus from the other diagrams as it should be. 



7 



b 



b 




Figure 2: Ghost two-point function of Landau-gaugo Yang-Mills theory. Solid lines denote gluons, dotted ones 
denote ghosts. The gray box represents the regulator insertion and blobs are dressed vertices. 



2.3. Derivation of Dyson-Schwinger equations 

For the sake of completeness we sketch the derivation of Dyson-Schwinger equations; details 
can be found in Ref. |98| . For a short description of a graphical derivation we refer the reader to 
Ref. [40j. To this end, we start with the following integral of a total derivative which yields zero: 



0= J D[q^y 



-S+<pjJj 



= m 



ss 



J,. 



SS 



4,'.=S/5J, 




Z[J]. 



(21) 



We are allowed to interchange integration and differentiation when we replace the fields by deriva- 
tives with respect to the corresponding sources. This yields the DSEs for full correlation functions. 
Next, we replace Z[J] by e^^'^l and multiply the resulting equation with e^^^'^l from the left. 
Using 



-WIJ] 



5W[J] 
6J, 



S 

SJ,, 



(22) 



we find 



5S 



(23) 



_ S\V[J] 



S.Ji ~^ SJi 



This is the (functional) DSE for connected correlation functions. To obtain the DSE for IPI ftmc- 
tions we perform a Legendre transformation of W with respect to all sources. Using 5W[J]/ 5Ji — 
and 



we are led to 



5J^ 



SJi (5$, 



S SW S 
SJi SJj S^j 



S'^W s 
SJiSJj (5$j 



(24) 



ST 

s¥,. 



5S 



(25) 



From this equation it is straightforward to derive DSEs for all IPI functions by taking derivatives 
with respect to the fields. Recall that the summation over j includes a sum over all fields of the 



theory. The command doDSE evaluates Eq. ( 25 ) and computes the n-point functions. For explicit 
examples we refer the reader to Ref. [55] . 



3. Basic usage of DoFun 

3.1. Actions 

The package DoDSERGE represents the heart of DoFun. It contains the two commands doDSE 
and doRGE which allow the derivation of DSEs and RGEs, respectively. As input these commands 
use a symbolic form of the physical (effective) action. Contrary to the physical action written 
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down with all indices and momenta, we shall refer to this symbolic form as symbolic action in 
the following. The symbolic action is a list of the types of fields and their interactions. This 
already suffices to determine the structure of the equations in terms of generic rt-point functions 
and propagators. 

No specific information about the fields like their nature (e. g. scalar, spinor or vector field) 
appear in the definition of the symbolic action. The user only needs to provide the information 
whether a field is real, complex or Grassmann-valued. Details, such as the color or spin structure, 
only need to be provided when the package DoAE is used. In fact, the command getAE allows to 
generate specific expressions from the generic correlation functions obtained from doDSE and doRGE, 
respectively. The resulting expressions can then be simplified by, e. g., contracting indices. DoFun 
offers only a tool to perform the simplest kind of contractions, namely those involving Kronecker 
deltas. The implementation of an advanced handling of indices such as Lorentz or color indices is 
left to the user. However, we would like to stress that programs for these tasks exist and can be 
applied to the output of DoFun. 

While the starting point for both doDSE and doRGE is an action, there is an important difference: 
doDSE requires the microscopic action S*, whereas doRGE is based on an ansatz for the effective 
average action Ffe. The latter is related to the classical action by fc — > A: Vk^A — S. 

As a simple example we consider a ip'^ theory. The symbolic action for this theory reads 

actionS={{(^, (^}, {ip , ip , ip , ip}} 
This action represents the bare action 

I [ cf^q d'^r d'^s ^ , , ^ , , , , 
+ 4! I (2.).(2.).(2.).^^('^M^M^M-'^-^-^) (26) 

for doDSE with A being the coupling constant. For doRGE, the symbolic action actions represents 
the ansatz for the effective average action 

TfcM = j ^^(/.(g)Z,(92)gV(-9) (27) 

1 f d'^q d'^r d'^s 
^4!i (2^ {27:)d {2Tr)d ^"^'^^ r,s,-q-r- sMq)p{r)ip{sM-q - r - s) 

where Zk{q^) is the so-called wave-function renormalization and is a momentum-dependent 
coupling constant. Note that in the definition of the symbolic action actions no signs or numerical 
factors appear. The program computes the latter from the multiplicity of the fields and their 
statistics. The bosonic or fermionic nature in this simple example is determined by the two-point 
function: if the list entries of length two have twice the same entry, e.g. {(^, p)} in the example 
above, the field is considered to be bosonic. If there are two different fields, e. g. {c, c}, then they 
represent as Grassmann-valued fields where c is the anti-field of c. 

There are two cases when we have to override this rule: (a) the fields mix at the two-point 
level as it is the case, for instance, in the Gribov-Zwanziger formalism |117j . (b) we have complex 
bosonic fields. This can be done with the option specif icFieldDefinitions, see below for a 
specific example. 

We would like to stress that the action actions can also stand for a different theory, for 
example, 

AY.j§f.i!^^.^MWWi>)9>(^,-r-sy (28) 

This ambiguity of the symbolic action is due to the fact that the combinatorics of the equations 
can be the same for various theories. 

As we use an ansatz for the effective average action in the derivation of RGBs, we may include 
operators which are not present in the classical action S. These operators may be generated in the 
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RG flow due to quantum fluctuations. In any case, the initial conditions for the RG flow equations 
are chosen such that F/j — >■ S* for fc — > A. Thus, the ansatz for the effective average action together 
with the initial conditions determine the quantum field theory under consideration. For instance, 
we may use the following symbolic action for a one-component scalar field theory which is invariant 
under the transformation — —ip: 

However, we may also use 

{{(^, if}, {^p, if, p, p}, {ip>, p, p, p, p, p}} 

The initial value of the six-boson interaction is then set to zero at the UV scale fc = A to ensure 
F/j — >■ iS* for A: — >■ A. Alternatively, the action under consideration can be specified by the maximal 
order in a given type of field, e. g. 

and 

6}} 

respectively. To check this we use the command generateAction which generates the correspond- 
ing action. For 

generateAction [{{(^, 4}}] 

we find that the action consists of three terms, namely a two-, three- and four-point function. We 
can enforce the reflection symmetry ^ — by adding the argument even: 

generateAction 4, even}}] 

Only vertices with an even number of legs are now included in the action. 

Let us now discuss complex scalar flelds. To define the symbolic action, we have to use the 
option specif icFieldDefinitions. Using 

actionCS={{(^, .^}, {p,p,p,p}} 

the commands doDSE and doRGE assume that p is a Grassmann-valued field by default. To avoid 
this, we use the option specif icFieldDefinitions in the derivation of functional equations. For 
the present case it reads 

specif icFieldDefinitions -> {p, ip} 

and is given as an argument to doDSE or doRGE. An example is provided in Section |4.1[ 

For theories in which fields mix at the level of the two-point function, such as the local Gribov- 
Zwanziger action [11 7| , we refer to a tutorial provided in the Documentation Center of Mathemat- 
ica. Indeed DoFun has already been applied successfully to such types of theories [HSl 1 1041 1106) . 

In order to avoid confusion which types of ansatze for the effective average action, viz. which 
expansions, can be used with DoFun we want to state explicitly that all ansatze are suitable as long 
as n-point functions can be defined from them. This includes also the usual derivative expansion. 

3.2. The basic constructs of DoFun 

DoFun uses the command generateAction to translate symbolic actions given in list form into 
its 'internal' language. The basic quantity of this language is the op construct. It acts as a container 
for fields, propagators, regulator insertions and vertices and as such it is responsible for keeping 
order in the expressions, op constructs also represent the summation/integration over indices 
and they are indispensable for the implementation of anti-commuting fields. The most general 
expression appearing in DoFun consists of a sum of op functions with an arbitrary complex-valued 
number as coefficienip^ Any numerical factors inside an op operator are immediately put in front of 
it. Before we discuss more specific examples, we introduce fields, propagators, regulator insertions 
and vertices. 



-'^"Notc that fields, vertices, regulator insertions or propagators are not permitted as coefliicients but only inside 
op functions. 
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Figure 3: Two (simple) examples for loop diagrams as they appear in Dyson-Schwinger equations, sec text for 
details. Solid lines represent a propagator, a blob a dressed vertex and the circle an external field. 



Quantity 


Meaning 


{A,i} 


field A with index i 


P[{A,i}, {A,j}] 


propagator of field A from i to j 


V[{A,i}, .... {B,j}] 


n-point function with legs {A,i}, . . . , {B, j} 


S[{A,i} {B,j}] 


bare n-point function with legs {A,i}, . . . , {B, j} 


dR[{A,i>, {A,j}] 


regulator insertion 


op[. . .] 


contains fields, propagators, 
regulator insertions and vertices 



Table 1: Overview of basic quantities of DoDSERGE. 



• Fields are represented by a simple list, e.g. {A , i}, where the first entry is the label of the field 
and the second entry its generic index. At this level we only use one index which comprises 
all physical indices including momenta. Fields and indices represent the two smallest bits of 
information. 

• In DoFun propagators and vertices only have fields as arguments. However, both propa- 
gators and fields are allowed as arguments of the op construct. A propagator is defined as 
P[{A,a},{B,b}] where the lists {A, a} and {B,b} denote its fields. 

• In the context of functional RG equations a regulator insertion is represented by dR with two 
fields as argument, e. g. dR[{A,a},{B,b}] . 

• In the framework of DoFun we have to deal with bare and dressed n-point functions denoted 
by S and V, respectively. The arguments of S and V denote the external legs of the vertices. 
For RGBs only dressed n-point functions V appear. 

While for bosons the order of fields is irrelevant, it is crucial for fermions. We have summarized our 
conventions for fermions in [Appendix A[ With these ingredients it is possible to draw Feynman 
diagrams with DoFun. For example: 

op[S[{W,i>,{W,k},{W,l>] , P[{W,k},{W,ks>] , P[{W,l>,-CW,ls}] , 
V[{W,ks},{W,ls},{W,j}]] 

see left panel of Fig. [3] Here W is a real bosonic field and the external indices of the diagram are i 
and j . A tadpole diagram with an additional external leg is given by 

op[S[{W,i>,{W,k},{W,l>,{W,m},{W,j}] , P[{W,k},{W,l>] , -CW.m}] 

see right panel of Fig. [3j Here, the field argument {W,m} of the op function denotes an external 
field. A summary of the basic objects of DoFun can be found in Table [T] 



11 



Vc/ LLcXiiLiiLV 


N/l o PI n m o' 


A [mom , i , j , . . . ] 


field A with momentum mom 
and indices i, j , ... 


P [A [moml . . .] ,A [mom2 ,i2,j2, . . .]] 


propagator oi neld A 

with indices il, i2, j 1, j2, ... 


V [A [moml ,il,...],...,B [momn , in , . . . ] ] 


71- point function with legs 
of nefds A, . . . , B and 
indices il, . . . , in 


S LA Lmoml ,il,...J,...,B Lmomn , in , . . . J J 


bare n-point function with legs 
of fields A, . . . , B and 
indices il, . . . , in 


dR[A[moml,il, jl, . . .] ,A[mom2,i2, j2, . . .]] 


regulator insertion 


op[. . .] 


contains fields only 



Table 2; Overview of basic quantities of DoAE and DoFR. 



3.3. Feynman rules 

In order to obtain physically meaningful expressions from the output of doDSE and doRGE we 
need to specify the actual (physical) meaning of the objects V, P, dR and S. Eventually this boils 
down to a specification of the Feynman rules. The definition of these objects is left to the user, i. e. 
the user has to provide the details concerning the propagators and the vertices. The Feynman rules 
can either be derived by hand or with the aid of independent packages, such as FeynRules |118j . 
Depending on the package used to derive the Feynman rules, the user may need to adapt the sign 
convention in these rules to ensure that they are compatible with our conventions for the definitions 



of n-point functions, see Eq. (5b) 



Alternatively, it is possible to use the package DoFR, which is part of DoFun. This package 
also allows to derive Feynman rules from a given action. Independent of the method we choose, we 
need to define the propagators and vertices to further process the output of the commands doDSE 
and doRGE. In general, this is done by adding appropriate definitions of V, P, dR and S. To be more 
specific, a propagator in symbolic form is given by 

P[{A,i},{A,j}] 

Here, A can be any bosonic field. Let us now assume that A is a gluon field, i. e. it has Lorentz 
and color indices. The explicit representation of the propagator is then given by 

P[A[pl, mu, a], A[p2, nu, b] ] 

The arguments are the fields whose arguments are the momenta and the (color) indices. Here, mu 
and nu are Lorentz indices and a and b are adjoint color indices. In order to define the propagator 
we evaluate the following expression: 

P[A[pl_, mu_, a_] , A[p2_, nu_, b_] , explicit->True] : = 
delta[color , a,b] (metric [mu,nu] -pi [mu] pi [nu] /pl~2)/pl~2/Z [pl~2] 

At this point several comments are in order: 

• The option explicit->True is required to force DoFun to use this definition of the propa- 
gator, see also below the function getAE. 

• The representation of objects like Kronecker deltas, four-momenta or dressing functions is 
completely left to the user. For illustration we have used a self-evident notation here. For 
example, metric [mu,nu] denotes the metric tensor and pi [mu] is a four-vector. 

• The momentum conservation of the propagator is taken into account implicitly, i. e. we 
define the propagator as a function depending on two different momenta but we only use 
one of these arguments on the right-hand side. Note that the output of getAE already takes 
momentum conservation into account. 

Other quantities, such as the vertices V and S or the regulator insertions dR, are defined ac- 
cordingly, see table [2j For example, the dressed ghost-gluon vertex can be defined as follows: 
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V[A[pl_, mu_, a_] , cb[p2_, b_] , c[p3_, c_] , explicit -> True]:= 
1 g structureConstant [a,b, c] 
(ghgDressingl[pl,p2,p3] p2 [mu] +ghgDressing2 [pi ,p2 ,p3] pl[mu]) 

Next, we turn to the package DoFR which can be used to derive the Feynman rules. The use 
of the commands of this package requires a physical action expressed in terms of the op operator. 
Note that in this case the syntax of the op operator differs from the syntax which is used to define 
a symbolic action: Only fields are permitted here as an argument for the op operator. Moreover 
they have to be given in an explicit form: 

field [momentum, indexl, index2, ...] 

where field is the name of the field. Its first argument is always the momentum and the others 
correspond to the labels of the indices of the fields, e. g. color or Dirac indices. The order of 
the indices is determined by the user. Note that the physical action always has to be given in 
momentum space. In general it is automatically assumed that there is an integral over reappearing 
momentum labels and a sum over reappearing indices. Integrals and sums are never written out 
explicitly in DoFun. As an example we consider the action of an iV-component scalar theory up 
to fourth order. First, we need to define phi as a field. In addition, we have to specify its indices. 



This is done with the command def ineFieldsSpecif ic, see Section 4.3 for more details: 
defineFieldsSpecific [{phi [momentum, type]}] ; 

DoFun now knows that phi is a bosonic field with one index type. We can check our specifications 
with 

{bosonQOphi, HeadOphi} 

— > {True, boson} 

and 

indices [phi] 
— > {type} 

The action in momentum space, 

WfWY'Wf^^'' r, - r - sW{qW{rW{sW{-q - . - ,s), (29) 

is then given by 

action=q~2 op[phi[q,i] ,phi [-q, i] ] /2 
+lajiibda [q,r , s ,-q-r-s] /8 op [phi [q, i] , 
phi [r , i] , phi [s , j ] , phi [-q-r-s , j ] ] ; 

To define an action within DoFun we have to obey a few rules. For example, the employed 
dummy indices must be unique variables and the momenta must be recognizable as momenta. To 
facilitate this task we may use the command convertAction which attempts to rewrite an action 
given by the user into such a form: 

actionC=convertAction [action] 



— > 1/2 q$744~2 op [phi [q$744 , dummy [1]], phi[-q$744, dummy [1]]] + 
1/8 lambda [q$748 , q$749, q$750, -q$748 - q$749 - q$750] 
op [phi [q$748 , dummy [2] ] , phi [q$749 , dummy [2] ] , 
phi[q$750, dummy [3]], phi [-q$748 - q$749 - q$750, dummy [3]]] 

The employed momenta are uniquely labeled q$i. The internal indices in action have been 
replaced by the placeholders dummy [ . . . ] which represent unique dummy indices. 

Finally we can use the command getFR to derive the two-point and four-point functions: 
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twoPoint=getFR[actionC, {phi[pl, i] , phi[p2, j]}] 



— > pl"2 delta[type, i, j] deltam[pl + p2] 

Here, the S distribution in momentum space is represented by 
deltam[pl+p2] 

which corresponds to (27r)'*(5(pl+p2) and delta [type, i, j] denotes a standard Kronecker delta 
with indices i and j . We shall return to the function delta below. 
The general syntax of getFR is 

getFR[actionC, {fieldl [momenta and indices], 
fields2 [momenta and indices], . . . }] 

From this we infer how to get the four-point function: 

f ourPoint=getFR [act ionC , 
{phi [pi, i] , phi[p2, j] , plii[p3, k] , phi[p4, 1] >] ; 

In the present example the four-point functions consists of 24 terms. This is the most general 
case but usually an approximation suffices. For instance, for a study of an 0{N) model a point- 
like approximation to the quartic coupling may already be sufficient to capture a wide range of 
physics}^ In this case we are left with three terms: 

fourPoint/ .lambda [ ] :>lcmibda 

— > -Icunbda delta [type, i, 1] delta [type, j, k] 

deltam[pl + p2 + p3 + p4] - 
lambda delta [type, i, k] delta [type, j, 1] 

deltam[pl + p2 + p3 + p4] - 
lambda delta [type, i, j] delta [type, k, 1] 

deltam[pl + p2 + p3 + p4] 



The minus sign is due to our definition of the vertices, see Eq. (5b). 

Having derived the equations for the two- and four-point functions we now have to define the 
specific form of the propagator and the vertex functions P and V, respectively. We infer the form 
of the propagator from twoPoint, which has been calculated above: 

P [phi [pl_ , indl_] , phi [p2_,ind2_] , explicit->True] : = 
delta [type , indl , ind2] /pl"2/Z [pl~2] ; 

where we added the dressing function since P represents the full propagator. It is important to 
stress that no S distribution should appear in the definition of the propagator. The latter is 
implicitly assumed. Similar to the propagator we can define the vertex. While the propagator is 
defined to be the inverse of the two-point function, we can directly employ the output of getFR to 
define the vertex function: 

V[phi[pl_, i_] , phi[p2_, j_] , phi[p3_, k_] , phi[p4_, 1_] , 
explicit -> True] = 
getFR [actionC , 

{phi [pi, i] , phi[p2, j] , phi[p3, k] , phi[p4, 1]}]/ 
deltam[pl + p2 + p3 + p4] // Expand; 

We would like to stress again that no momentum-conserving S function should appear in the 
definition of the vertex functions. Bare vertices S and regulator insertions dR are defined in the 
same way as the vertex functions. 

Finally we would like to comment on a function offered by DoFun which may prove useful in 
many cases, namely delta. It is a generalization of the standard Kronecker delta which allows 
to relate the indices to a certain type of object. For example, delta [ind, i, j] is a Kronecker 
delta with indices i and j associated with the index type ind. Let us give a few basic examples 
for the application of the function delta: 



^In an RG approach the quartic couphng is a scale-dependent quantity, even in the point-like limit. 
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delta[l,l] — > 1 
delta [0,3] — > 
delta[ind,a,a] — > dim[ind] 

where dim represents the dimension of the representation associated with this index. Ahhough it is 
not mandatory to specify the type of object to which the indices belong, we highly recommend to 
do so in more involved cases with several different indices and types of objects. With the function 
integrateDeltas DoFun offers the possibility to 'integrate out' these delta functions. Usually 
this function even works in cases in which one of the indices is part of a different function, e. g. 

someFunction [a, b] delta[a, c] // integrateDeltas 
— > someFunction [c , b] 

Note that integrateDeltas checks whether an index appears more than twice in delta functions. 
The functions delta and integrateDeltas arc also used by getFR. In the subsequent sections we 
discuss the usage of the main functions of DoFun by means of specific examples. 

4. Derivation of functional equations 

Derivation of Dyson-Schwinger equations 
The function for the derivation of DSEs is doDSE. In the simplest case the function call reads 

doDSE [action, derivatives] 



The main input is the symbolic action which is called action in the example abovq^ see Section 3.1 
for details. The argument derivatives is a list of fields corresponding to the legs of the n-point 
function we would like to compute. A simple example is the DSE for the two-point function of a 
scalar theory, see also Fig. |4} 

doDSE [{{phi ,phi} , {phi , phi , phi , phi}} , {phi, phi}] 

--> op [S [{phi, il}, {phi, 12}]] - 

1/2 op [S [{phi, 11} , {phi, 12}, {phi, rl}, {phi, si}], 

P[{phi, rl}, {phi, si}]] - 

1/6 op [S [{phi, 11}, {phi, rl}, {phi, r2}, {phi, si}], 

P[{phi, rl}, {phi, s2}] , P[{phi, r2}, {phi, t2}] , 

P[{phi, si}, {phi, u2}] , 

V[{phi, 12}, {phi, s2}, {phi, t2}, {phi, u2}]] 

When we do not attach indices to the fields in the list derivatives, the indices are attached auto- 
matically to the legs by doDSE. In the present example, doDSE has chosen 11 and 12 as index labels. 
Alternatively, it is possible to choose specific indices by including them in the list derivatives: 

twoPointPhi=doDSE [{{phi ,phi} , {phi , phi , phi , phi}} , {{phi, a}, {phi,b}}] 

— > op [S [{phi, a}, {phi, b}] ] - 
1/2 op[S[{phi, a}, {phi, b}, {phi, rl}, {phi, si}], 

P[{phi, rl}, {phi, si}]] - 
1/6 op [S [{phi, a}, {phi, rl}, {phi, r2}, {phi, si}], 

P[{phi, rl}, {phi, s2}] , P[{phi, r2}, {phi, t2}] , 

P[{phi, si}, {phi, u2}], 

V[{phi, b}, {phi, s2}, {phi, t2}, {phi, u2}]] 

For a fermionic theory the derivation of functional equations works in the same way. Consider 



^^In principle, doDSE and doRGE also accept physical actions. However, the physical actions arc transformed into 
symbolic actions by doDSE and doRGE . 
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Figure 4: The two-point function DSE of 95'* theory. 




Figure 5: The two-point function DSE of (V"/")^ theory. 



twoPointPsi=doDSE [{{psi , psib}, {psib, psib, psi, psi}}, 
{{psib, a>, {psi, b»] 

— > op [S [{psi, b}, {psib, a}]] + 
op[S[{psib, a}, {psib, rl}, {psi, si}, {psi, b}] , 

P[{psi, si}, {psib, rl}]] - 
1/2 op[S[{psib, a}, {psib, rl}, {psi, r2}, {psi, si}], 
P[{psi, s2}, {psib, rl}] , P[{psi, r2}, {psib, t2}] , 
P[{psi, si}, {psib, u2}] , 

V[{psib, t2}, {psib, u2}, {psi, s2}, {psi, b}] ] 

Here, psi is the fermion and psib the anti-fermion. The first entry in the list action corre- 
sponds to the kinetic term whereas the second term corresponds to a four-fermion interaction, 
see also Sect. |5.2| The result is depicted in terms of Feynman diagrams in Fig. [Sj The notation 
for fermions as {fermion, anti-fermion} is quite convenient and is used throughout DoFun. 
However, the treatment of complex bosonic fields also requires different labels for the fields. To 
avoid the identification of complex bosonic fields as Grassmann fields we have to use the op- 
tion specif icFieldDefinitions in doDSE. To be specific, we consider an action which describes 
a toy model of a complex scalar field phi and a ferniionic field psi: 

twoPointMixed=doDSE [{{phi , phib}, {psi, psib}, 
{psib, psi, phib, phi}, {phib, phib, phi, phi}, 
{psib, psib, psi, psi}}, 
{{phi, b}, {phib, a}}, 

specif icFieldDefinitions -> {phi, phib, {psi, psib}}]; 

The arguments of specif icFieldDefinitions are the bosonic fields and a list of the fermion pair. 
The output is shown graphically in Fig. [6] 

^.2. Derivation of renormalization group equations 

The derivation of RGEs with the command doRGE works very similar to the derivation of DSEs. 
We start with two simple examples: the fiow equations for the effective average action itself and 
the two-point function for a iy9^-theory, see Fig. [Tj 

doRGE[{{phi, phi}, {phi, phi, phi, phi}}, {}] 

--> 1/2 op[dR[{phi, rl}, {phi, si}], P[{phi, si}, {phi, rl}]] 

and 



-1 b ^ -1 _ ^ f \ b 



b a . b 



Figure 6: The two-point function DSE of a theory with complex scalars ip and fermions i/i with interactions (ipip)^, 
i'tpip)'^ and (p^p'ipip. Bosons are denoted by a continuous line, fermions by a dotted one. 
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Figure 7: The flow equations for the effective action (left) and the two-point function (right). 



doRGE [{{phi , phi}, {phi, phi, phi, phi}}, {phi, phi}] 

--> 1/2 op[dR[{phi, rl}, {phi, si}], P[{phi, tl}, {phi, rl}], 
P[{phi, si}, {phi, vl}] , 

V[{phi, i2}, {phi, il}, {phi, vl}, {phi, tl}]] 

The rules for fermions and complex fields are the same as for the derivation of DSEs. 

For doRGE we have the additional option tDerivative. It allows to suppress the derivative 
dt which attaches the derivative of the regulator to the internal lines of the IPI diagramj^ 
This is useful if one is only interested in the structure of the RG equations. For illustration, 
we show the output of doRGE for the three-point function of a (p^ theory in Fig. [s] as obtained 
for tDerivative->True (top panel) and for tDerivative->False (bottom panel). To obtain the 
output shown in the top panel of Fig. |8] we have used 

threeR = doRGE[{{phi, phi}, {phi, phi, phi}}, {phi, phi, phi}]; 

whereas we have employed 

threeRnot = doRGE[{{phi, phi}, {phi, phi, phi}}, {phi, phi, phi}, 
tDerivative -> False] 

--> op [V [{phi, il}, {phi, rl}, {phi, r2}], P [{phi , r2}, {phi, s2}] , 
V[{phi, 12}, {phi, s2}, {phi, t2}] , P[{phi, t2}, {phi, u2}] , 
V[{phi, 13}, {phi, u2}, {phi, v2}] , P[{phi, rl}, {phi, v2}]] 

to obtain the output shown in the bottom panel. Here, we have not given the non-graphical output 
of doRGE for the case with tDerivative->True since it is rather lengthy. Note the different signs 



appearing in the two graphical representations which arise due to the derivative Eq. (151 



4.3. (Graphical) Representation of the output of doRGE and doDSE 

The output of doDSE and doRGE in terms of op functions does not look very transparent from 
a physical point of view. In the previous sections we have already made use of the fact that the 
output of doDSE and doRGE can be illustrated in terms of Feynman diagrams. Before we discuss 
this option in more detail, we would like to discuss a further possibility to represent the output of 
doDSE and doRGE within DoFun. 

The function shortExpression, or equivalently sE, transforms the fields, propagators, regulator 
insertions and vertices in the output of doDSE and doRGE into more familiar expressions. For 
example, a propagator is denoted by A. Consider the expression: 

ShortExpression [op [S[{W,i},{W,k},{W,l}] , P[{W,k},{W,ks}] , 
P[{W,l},{W,ls}] , V[{W,ks},{W,ls},{W,j}]]] 



This yields 



oirlsl Y-tlvlj Arltl \slvl 



^•^Note that diagrams with the same topology but the regulator insertion attached to different internal lines appear 
on the right-hand side of the functional RG equations. 
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Figure 8: Output of doRGE for the three-point function of a ip"^ theory in terms of Feynman diagrams as obtained 
for tDerivative->True (top panel) and for tDerivative->False (bottom panel). Note that the derivative is not 
indicated explicitly on the right-hand side of the bottom panel though it is contained in the expression. 



Symbol 


Standard value 


$propagatorSymbol 


A 


$bareVertexSymbol 


S 


$vertexSymbol 


r 


$regulat or Insert ionSymbol 


R 



Table 3: Symbols used in shortExpression. 



The field types are given as subscripts and their (collective) indices (momenta, . . . ) are given 
as superscripts. Propagators are represented as A and bare and dressed vertices as S and F, 
respectively. A regulator insertion is denoted by dtR- Note that these string objects can be 
changed with the aid of the variables given in table [Sj Moreover shortEexpression accepts style 
options, e. g. 

ShortExpression [op [S[{W,i},{W,k},{W,l>] , P[{W,k},{W,ks}] , 
dR[{W,ks},{W,ms}] , P[{W,ms},{W,m}] , P[{W,l>,-CW,ls}] , 
V[{W,m},{W,ls},{W, j}]] , FontSize->20, Bold] 

It yields 

fl TDsltl Qir l7-2 -pulvlj A rlsl Ar2vl a tlul 

'^t^ww'-'www^ www^ww^ww^ww 

Of course, the representation of functional equations in terms of Feynman diagrams is most 
intuitive from a physical point of view. To this end, we have included two commands, namely 
DSEPlot and RGEPlot. As these functions are very similar we only discuss RGEPlot here. Differ- 
ences exist in the output, e. g. RGEPlot adds a dt to the left-hand side of the functional equation. 
Usually RGEPlot is called as follows: 

RGEPlot [output of doRGE, style definitions for the fields] 

The first argument is the output of doRGE. Of course, modifications by the user are allowed as 
long as the basic structure, namely a sum of op functions, is not changed. The second argument 
is optional. It allows to determine the graphics style for the fields. It has to contain a list of lists: 
The first argument of each sublist gives the name of the field and the remaining arguments can 
be graphics directives specifying the style of the propagators of that field. For the graphics style 
of the anti-fermions RGEPlot automatically uses the same graphics style as for the fermions. Style 
definitions can be, for example, colors or directives like Dotted or Dashed. If we do not set the 
style of the fields explicitly, the propagators and vertices are labeled according to the names of 
their fields. Examples for these two cases are 

RGEPlot [ 
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a 




b 



Figure 9: Examples for the representations of RGEs in terms of Fcynman diagrams without (top) and with (bottom) 
user-defined styles for the fields. 



doRGE [ 

{{phi, phi}, {phi, phi, phi}, {phi, phi, phi, phi}}, 
{{phi, a}, {phi, b}}] 

] 

and 

RGEPlot [ 
doRGE [ 

{{phi, phi}, {phi, phi, phi}, {phi, phi, phi, phi}}, 
{{phi, a}, {phi, b}}] , 
{{phi, Black}} 

] 

The output is shown in Fig. [9] 

There are also options which allow to set the style of indices (indexStyle) or numerical coeffi- 
cients (f actorStyle) and to determine the number of diagrams shown in one row. These options 
are explained in detail in a dedicated part of the Documentation Center of Mathematica. Here, we 
only discuss the option output. It allows three settings: 

• List: Gives a list of individual graphs. 

• f orceEquation: Draws also the left-hand side of the equation. 

• complete (default): Gives the complete equation including the left-hand side if the expression 
contains several graphs. If the expression consists of a single graph, only this graph is shown. 

We would like to point out that we have to define the nature of the fields in order to use DSEPlot 
and RGEPlot. As discussed above, the fields are automatically defined when we use doDSE and 
doRGE but they can also be defined by hand with the aid of the command def inePields. Its 
syntax is 

def inePields [list of bosons, list of fermions, 
list of complex fields] 

The lists of fermions and complex fields are given in the usual double notation of fermions and 
anti-fermions, e.g., {{c,cb}, {q,qb}}. If there are no fields of a certain type, an empty list, 
{}, is required as argument. Note that the nature of the fields is defined automatically when 
generateAction is used. 
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Both functions, DSEPlot and RGEPlot, are based on the Mathematica function GraphPlot. 
This has the advantage that we can directly use a built-in function of Mathematica. On the 
other hand, there are a few drawbacks. For example, we are not aware of a simple way to plot 
wiggly lines to, e. g., represent gluons. Also, for higher vertex functions it might be the case 
that different lines appear on top of each other. Thus, the output looks wrong at first glance. 
Non-planar diagrams also represent a problem for DSEPlot, see a corresponding notebook in the 
Documentation Center. In general, however, DSEPlot and RGEPlot provide useful output which 
allows to check the associated functional equations. 

4.4- Algebraic expressions 

As illustrative the graphical representation may be, the full analytic expressions containing 
all indices and momenta are required for an evaluation of the functional equations. To this end, 
we have included the package DoAE. To translate the output of doDSE or doRGE into algebraic 
expressions we use the function getAlgebraicExpression, or short getAE. It uses the symbolic 
output of doDSE or doRGE and adds the indices and momenta as indicated by the user. Its syntax 
is 

getAE [exp, external momenta and indices, options] 

The first argument is an expression in symbolic notation, not necessarily the output of doDSE or 
doRGE. If it is a sum of several op-functions, getAE returns a list so that terms can be traced back 
to their origin. To obtain the sum instead, one uses the command Plus@@get AE [ . . . ] . The second 
argument is a list of the external legs together with labels for their momenta and indices. For 
every external leg there has to be a list entry as follows 

{field, generic index, momentum, real indices} 

Here, generic index refers to the index label which appears in the output of doDSE or doRGE. 
Recall that the indices are chosen according to il, 12 and so on, if we do not explicitly provide 
index labels, momentum can be a symbol, e. g., pi, or a number, real indices are the indices 
corresponding to this leg. Note that Fcynman rules have to be properly defined by overloading the 



functions P, V, S and dR as described in Section 3.3 



The command getAE replaces the generic indices appearing in the output of doDSE and doRGE 
with conventional symbols for the momenta and indices of the fields and applies the Feynman 
rules for the propagators and vertices. Therefore we need to define the indices with the function 
def ineFieldsSpecif ic: 

defineFieldsSpecific [{boson [momentum, indexl, ...], ... 
{fermion [momentum, index2, ...], 
Einti-f ermion[momentum, index2, ...]}, ...] 

The argument is a list of fields where bosons and fermions are specified in the typical manner of 
grouping fermion and anti-fermion into sublists. The arguments of the fields denote the labels 
of their indices, momentum is the obligatory first argument followed by the other indices. The 
indices of a field can be checked with the command indices. Here we do not need to worry about 
complex fields. Note that defineFieldsSpecific should not be confused with defineFields 
which is required for plotting graphs. 

The names of dummy indices depend on the specific type of index. For example, Lorentz 
indices, labeled lor, are named \ [Mu] , \ [Nu] and so on by getA^^ Apart from Lorentz indices 
color indices in the adjoint representation, labeled adj , are the only predefined indices. If not 
specified otherwise, the names a, b, c are assigned automatically to new index labels. Alternatively, 
it is possible to assign labels to specific index types with the aid of the function addlndices. In 
fact, this function assigns a list of labels to a specific index name, e. g. flavor: 

addlndices [{flavor , {i, j, k, 1, m, n}}] 

With resetlndices [] we can reset the index labels to the default definitions of DoFun. 

We illustrate the use of getAE with the aid of a comparatively simple theory, namely an A^- 
component scalar theory truncated at the level of the four-point function. Many of the required 



Note that the internal routine ensures that these names are unique. 
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steps to derive the functional equations of this theory have already been explained above. Therefore 
we only give the commands to derive the flow equations of the two- and four-point functions. First 
we define the fields 

defineFieldsSpecific [{phi [momentum, type]}] ; 

Next we specify the name of the indices. We choose i, j and so on: 

addlndices [{type, {i, j, 1, m, n>}] ; 

The Feynman rules for the propagator, the regulator insertion and the quartic vertex with a 
momentum-independent coupling lambda are defined as follows: 

P[ph.i[pl_, i_] , phi[p2_, j_] , explicit -> True] : = 

delta[type, i, j]/(Z[k, pl~2] pl~2 + R[k, pl"2]) 
dR[phi[pl_, i_] , phi[p2_, j_] , explicit -> True] := 

delta[type, i, j] dR[k, pl~2] 
V[phi[pl_, i_] , phi[p2_, j_] , phi[p3_, k_] , phi[p4_, 1_] , 
explicit -> True] := 
-Icimbda delta [type, i, 1] delta [type, j , k] - 
lambda delta [type, i, k] delta [type, j , 1] - 
lambda delta [type, i, j] delta [type, k, 1] 

To derive the RG equations for the two-point and the four-point functions we use 

twoR = doRGE[{{phi, phi}, {phi, phi, phi, phi}}, 

{{phi, a}, {phi, b}}] ; 
fourR = doRGE[{{phi, phi}, {phi, phi, phi, phi}}, 

{{phi , a} , {phi , b} , {phi , c} , {phi , d}}] ; 

Now we employ getAE and choose the external momenta to be zero: 

Plus(§(SgetAE[twoR, {{phi, a, 0, i}, {phi, b, 0, j}}] ; 
PlusOOgetAE [fourR, {{phi, a, 0, i}, {phi, b, 0, j}, 
{phi, c, 0, 1}, {phi, d, 0, m}}] ; 

We applied Plus to get the sums of the expressions. To bring the results into a convenient form 
we multiply them with appropriate Kronecker J's and sum over the corresponding indices. For the 
two-point function we obtain 

integrateDeltas [ 

delta [type, j, 1] delta [type, j, 1] Plus OQ 

getAE[twoR, {{phi, a, 0, i}, {phi, b, 0, j}}]] /. 
dim [type] : > Ntype 

~> -((lambda delta[type, i, j] dR[k, ql~2])/ 
(R[k, ql"2] + ql-2 Z[k, ql~2])-2) - 
(lambda Ntype delta[type, i, j] dR[k, ql~2])/ 
(2 (R[k, ql-2] + ql-2 Z[k, ql-2])-2) 

For the four-point function we find 

integrateDeltas [ 
delta [type, i, 1] delta [type, j, 1] delta [type, 1, 1] 
delta[type, m, 1] PlusOOgetAE [fourR, 
{{phi, a, 0, i}, {phi, b, 0, j}, {phi, c, 0, 1}, 
{phi, d, 0, m}}] 
] /. dim [type] :> Ntype 

— > -((24 lambda"2 dR[k, ql-2])/(R[k, ql"2] + ql"2 Z[k, ql"2])-3) - 
(3 lambda"2 Ntype dR[k, ql-2])/(R[k, ql-2] + ql-2 Z[k, ql-2])-3 

Specifying the regulator R[k,ql~2] and its insertion dR[k, ql~2] and adding the integral yields 
then the final results suitable for an analytic or numerical evaluation. 

This was a first quick introduction into the usage of getAE. For a more detailed introduction we 
refer the reader to the subsequent section. There, we also discuss how the equations for a regime 
with spontaneously broken 0{N) symmetry can be derived. 
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5. Examples 

5.1. 0(N) models 

0{N) symmetric scalar field theories play a very prominent role in theoretical physics. First 
of all, they represent a valuable testing ground for a study of spontaneous symmetry breaking. 
In particular, the critical exponents at the phase transition are accessible to a large variety of 
methods. Therefore studies of the critical behavior of 0{N) models allow us to benchmark the- 
oretical approaches. Within the functional renormalization group approach critical behavior has 
been indeed studied in great detail, see e. g. Refs. 

In this section we use several symbolic actions representing different levels of truncations. First, 
we define an ansatz for an action which contains interactions up to the ten-boson coupling. For 
the flows calculated below this effectively corresponds to no truncation: 

actionONSymbolic = {{phi, 10, even}}; 

We then specify a truncated action which only includes a quartic interaction term: 
actionONSTrunc2 = {{phi, 4, even}}; 
The fields are specified with one index type: 
defineFieldsSpecific [{phi [momentum, type]}] ; 

We define the corresponding physical actions below. However, we would like to add that it is 
possible to use a potential U{(j3^) to define the action instead of a sum of various interaction terms, 
e. g. U{(l)^) — (l/2)m^0^ -I- (A/8)0'* where t/)'^ = , . . . , ^at). This can be done as follows: 

actionONPotential = convertAction [ 
1/2 (q-2 Z[q-2] + R[k]) op [phi [q, a], phi[-q, a]] + U[phi]]; 

DoFun recognizes U [phi] as a self-interaction potential and its derivatives are denoted by der [ 
U[phi], {phi[pl,il], ...}]. 

For the representation of the internal lines in Feynman diagrams we choose a solid black line: 

fieldStyleON = {{phi, Black}}; 

5.1.1. Symmetric regime 

In case of unbroken 0{N) symmetry in the ground state we have (0) — 0. We may then use 
the following ansatz for the effective potential and the wave-function renormalization 

OO T 

Uk{p) = C/fe(0) + J2 ^Z'"' (30) 
^ — ' ml 

m—l 
OO „(m) 

Zkip) - E ^z^"' (31) 
^ — ' m! 

where p = 1/20^. 

The symbolic forms of the flow equations for Uk{0), the two-, four- and six-point functions are 
obtained as follows: 

zeroPoint = doRGE [actionONSymbolic , {}] ; 

twoPoint = doRGE [actionONSymbolic , {phi, phi}]; 

fourPoint = doRGE [actionONSymbolic , {phi, phi, phi, phi}]; 

sixPoint = doRGE [actionONSymbolic , {phi, phi, phi, phi, phi, phi}]; 



The results for the zero-, two- and four-point functions are depicted in Fig. 10 At this level of 
the derivation these equations are still exact. Higher rt-point functions can also be obtained in 
this way. However, the computing times increase due to the high number of terms arising at 
the intermediate steps of the derivation. For example, it takes roughly two minutes on an AMD 
Phenom II X4 quad core processor to derive the eight-point function; approximately 100.000 terms 
are generated during the derivation. At the end we are left with only 1.954 distinct diagrams. 
As a non-trivial check we have also derived the flow equation for the nine-point function. This 
involves about one million intermediate terms. As expected due to the underlying symmetries of 
the theory, these terms sum up to zero. 

Let us now define the propagator of the scalar field and the regulator insertion: 
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Figure 10: Flow equations for the effective average action and the two- and four-point functions. 



P[phi[pl_, a_] , phi[p2_, b_] , explicit -> True] := 

(delta [type, a, b])/(Z[pl~2] + R[pl~2, k] ) ; 
dR[phi[pl_, a_] , phi[p2_, b_] , explicit -> True] := 
delta[type, a, b] dR[pl~2, k] ; 

In the following we truncate the ansatz for the effective average action at the lowest non-trivial 
order, i. e. we only allow for a quartic interaction term. Thus the physical action reads 

actionOM2 = 

convertAction [ 
1/2 (p-2 Z[p-2] + R[k]) op[phi[p, i] , phi[-p, i]] + 
laiiibda2/8 op [phi [p, i] , phi[q, i] , phi[r, j] , 
phi[-p - q - r, j]]] ; 

The Feynman rule for the four-point function is then defined as follows: 

V[phi[pl_, i_] , phi[p2_, j_] , phi[p3_, 1_] , phi[p4_, m_] , 
explicit -> True] = getFR[actionON2, {phi [pi , i] , phi[p2, j] , 
phi[p3, 1], phi[p4, m] }] /deltam[pl + p2 + p3 + p4] // Expand; 

To obtain the algebraic expressions for the equations we use getAE. We start with the fiow 
equation for the effective average action: 

zeroPointAlg = Plus <§@ getAE [zeroPoint, {}] // integrateDeltas 

— > (dim[type] dR[ql"2, k])/(2 (R[ql'2, k] + Z[ql"2])) 

To derive the algebraic expressions for the flow equations of the two- and four-point functions with 
getAE we need to provide lists for the indices and momenta associated with the external legs: 

twoPointAlg = Plus @® getAE [twoPoint, 

{-[phi, il, pi, i>, {phi, 12, p2, j}>] // integrateDeltas 

~> -((lambda2 delta [type, i, j] dR[(pl + p2 + ql)"2, 

k])/((R[(-pl - p2 - ql)-2, k] + 
Z[(-pl - p2 - ql)-2]) (R[ql-2, k] + Z[ql-2]))) - 
(lambda2 delta[type, i, j] diin[type] dR[(pl + p2 + ql)"2, k])/ 

(2 (R[(-pl - p2 - ql)'-2, k] + Z[(-pl - p2 - ql)-2]) (R[ql~2, k] + 
Z[ql-2])) 
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and 



fourPointAlg = Plus 0® getAE [f ourPoint , 

{{phi, il, 0, i}, {phi, i2, 0, j}, {phi, i3, 0, 1}, 

{phi, i4, 0, m»] /. V[a___] :> /; Leiigth@{a} > 4 // 
integrateDeltas// Simplify 

— > -((lambda2~2 (delta [type , i, m] delta [type, j , 1] + 
delta [type, i, 1] delta [type, j , m] + 

delta [type, i, j] delta [type, 1, m]) (8 + dim [type]) dR[ql~2, 
k])/(R[ql-2, k] + Z[ql-2])-3) 

Note that we consider the point-Hke limit here, i. e. we set the external momenta to zero. Ac- 
cording to our truncation we also set the six-point function to zero with the replacement rule 
V [a ] : > / ; Length@{a} > 4. Alternatively, we could also define this vertex as 

V[phi[_], phi [_] , phi [_] , phi [_] , phi [_] , phi [_] , explicit->True] : =0 

Another way would have been to use the truncated action actionONSTrunc2 for the calculation 
of f ourPoint. The resulting expressions can be evaluated further. This is discussed for the case 
with broken 0{N) symmetry in the following section. 

5.1.2. (Spontaneously) Broken 0{N) symmetry in the ground state 

In case of a (spontaneously) broken 0{N) symmetry in the ground state the expectation value 
of the field is non-vanishing, (0) ^ 0, and the modes of the theory acquire different masses. In 
fact, we encounter one massive mode and A^— 1 Goldstone modes. For convenience, we expand the 
theory about {(j)) = (j)oS'^^ where we consider (po to be space-time independent. For simplicity, we 
shall study a truncation in which we only take into account the leading order terms in an expansion 
of the potential and the wave-function renormalization in powers of p = (l/2)(/)^ : 



At finite temperature the constant term in the expansion of the potential Uk is related to thermo- 
dynamic observables. The subscript denotes the vacuum expectation value, i.e., po :— {l/2)(f)Q. 
The corresponding physical action reads 

defineFieldsSpecific [{phi [momentum, type], phiO [momentum, type]}]; 
actionON2 = 
convert Action [ 
1/2 (p-2 Z[p-2] + R[k]) op [phi [p, i] , phi[-p, i]] + 
lambda2b/8 op[op[phi[p, i] , phi[q, i]] - 
op[phiO[p, i] , phiO[q, i]], 
op [phi [r, j] , phi[-p - q - r, j]] - 
op[phiO[r, j], phiO[-p - q - r, j]]]]; 

Note that we have defined phiO as a field. This allows us to use it in the derivation of the n-point 
functions. We have used the op function as multiplication operator and therefore op functions 
appear as arguments of an op functionFj The Feynman rules are given by 




P[phi[pl_, a_] , phi[p2_, b_] , explicit -> True] := 
delta [type, a, 1] delta [type, b, 1]/ 

(Z[pl~2] pl~2 + R[pl~2, k] + 2 rhoO lambda2b) + 
(delta[type, a, b] - delta[type, a, 1] delta[type, b, 1])/ 
(Z[pl-2] pl-2 + R[pr2, k]) ; 

dR[phi[pl_, a_] , phi[p2_, b_] , explicit -> True] := 



l^Note that op[. . . , op[a], ...]= op[..., a, ...]. 



C/fc(p) = f/fe(po) + y(p-Po)', 



(32) 
(33) 
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delta[type, a, b] dR[pl"2, k] ; 



V[phi[pl_, i_] , phi[p2_, j_] , phi[p3_, k_] , phi[p4_, 1_] , 
explicit -> True] = 

getFR[actionON2, {phi[pl, i] , phi[p2, j] , phi[p3, k] , phi[p4, 1]}, 

symmetry -> broken] /deltam [pi + p2 + p3 + p4] // Expand; 

We split the propagator into N — 1 massless modes (Goldstone modes) and one massive mode. 
For the sake of simphcity, we consider the wave-function renormahzation Z to be identical for the 
massive and the radial mode. In the derivation of the four-point function we used the option 
symmetry -> broken to indicate that the vacuum expectation value of the field does not vanish. 

With these Feynman rules we can employ the symbolic expression for the zero-point function 
from the previous section to obtain the flow equation for Uk{po)- 

zeroPointAlg = Plus OS getAE [zeroPoint , {}] // integrateDeltas 

--> -(dR[ql-2, k]/(2 (R[ql-2, k] + ql-2 Z[ql-2]))) + 
(dim[type] dR[ql'2, k])/(2 (R[ql~2, k] + ql"2 Z[ql'2])) + 
dR[ql"2, k]/(2 (2 lambda2b rhoO + R[ql~2, k] + ql'2 Z[ql~2])) 



The left-hand side of the equation reads 



dtTk = / d'^xdtUki^o) = dtUkiM / d'^x. (34) 

0=00 J J 

The integral over d'^x corresponds to the momentum conserving S distribution on the right-hand 
side, (27r)''5(0). The latter is not given explicitly in the output of doRGE but always implicitly 
assumed. Thus, the flow equation for Uk{po) reads 

aTT( \ I f dS dtRkjq^) 

^Zkiq^)q^ + Rkiq') + 2\2P0^' ^ ' 

This equation agrees with the flow equation found in Refs. [T^ [^ . 

For the derivation of the flow equations of po a-nd the quartic coupling A2 we introduce the 
renormalized dimensionless mass and coupling: 

K = Zkk^^'^pa , 
A2 = Z-^k''-^\2 . (36) 

Here, we consider a truncation in which we only include the flow of pq and the quartic interac- 
tion. Higher-order interactions, such as a six-boson interaction, are set to zero. Within DoFun it 
is convenient to start with the derivation of the flow equation of the highest-order coupling. This 
is due to the fact that the flow equation of the n-boson coupling depends directly on the flow of 
the (n -I- 2)-boson coupling via a term ^ dt(f>Q, see Ref. [22] and the discussion of the flow of po 
below for details. 

The left-hand side of the flow equation of the four-point function reads: 



d(l)'{x)S(j)^{y)S(j)''{z)(t)"^(z)^ 

= Zlk'^-'^iS'^S^"' + (5*'(5J™ + (5™,5^')(4 - d - 2?7 + dt)\26{x - y)6{x - z)S{x - u), (37) 

where rj — — 9(lnZ. The Fourier transformation for vanishing external momenta yields 

6^ 



(5(/)*(O)(5(/)J(O)J0'(O)0™(O) 

Zlk*-'^{S'^^"' + S'^S^"" + (5™(5^')(4 - d - 2ry + dt)X2 I d'^x. (38) 
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On the right-hand side we obtain the algebraic form as foUowi 



f ourPointTrunc2 = doRGE [actionONSTrunc2 , {phi, phi, phi, phi}]; 
f ourPointT2Alg = 
Plus 00 getAE[ 

f ourPointTrunc2 , {{phi, il, 0, i}, {phi, i2, 0, j}, 
{phi, i3, 0, 1}, {phi, i4, 0, m}}] // integrateDeltas ; 

where we derived the flow of the four-point function from actionONSTrunc2 in the first line. In 
particular for a study of critical phenomena it is convenient to introduce dimensionless quantities, 



see Eq. (36). We define the following set of rules to rewrite the output accordingly: 



dimLessRulesON = {R[ql'2, k] :> ql~2 r[y], 
rhoO :> Z[ql"2]~(-1) kappa k'(d - 2), 
lambda2b :> Z[ql~2]'2 lambda2 k~(4 - d) , 

dR[ql~2, k] :> ql"2 dtr [y] , ql :> Sqrt [y ] k, dql :> k~2/2/ql}; 

Multiplying f ourPointT2Alg with S'^^S^^S''^S"^^ and applying the rules dimLessRulesQN we find 

anglelnt[d_] :=4v[d]; 
f ourPointT2Algy = 
Simplify [(cinglelnt [d] dql ql~(d - 1) integrateDeltas [ 

delta [type, i, 1] delta [type, j, 1] delta [type, 1, 1] 
delta[type, m, 1] f ourPointT2Alg] //. dimLessRulesON), 
d \ [Element] Integers] 

— > 6 k'(4 - d) lambda2~2 y~(d/2) 

dtr[y] v[d] Z[k'2 y] '4 (l/(y"3 (r[y] + Z[k"2 y])'3) - dim[type]/ 
(y-3 (r[y] + Z[k-2 y])-3) - 

9/(y r[y] + (2 kappa lambda2 + y) Z[k'2 y])~3) 

Here we have introduced the integral measure angelint [d] dql ql" (d - 1) and v[d] represents 
Vd — (7r''/^2''+^r((i/2))^^. On the left-hand side the projection yields a factor of 3: 

integrateDeltas [ 

delta [type, i, 1] delta [type, j, 1] delta [type, 1, 1] delta [type, m, 
1] (delta [type, i, j] delta [type, 1, m] + 
delta [type, i, 1] delta [type, j , m] + 

delta [type, i, m] delta [type, j, 1])] //. dimLessRulesON 

--> 3 

Combining both sides yields 

d,X,={d-^ + 2r,)X, + 2vAl r dy //^ ({N - 1) t^^^^ 

Jo \ (^fc + r(y))3 

Zldtr(y) 



(39) 



{Zk{y + 2KX-2) + yr{y))\ 
Recall that doRGE computes the n-point function which is defined to be the negative derivative of 



the effective average action for n > 2, see Eq. (5b) 



Next, we derive the algebraic form of the fiow equation for the two-point function in the limit 
of vanishing external momenta. To this end, we need to rederive the expression for the two-point 
function since now additional contributions due to the broken 0{N) symmetry appear on the 
right-hand side. The result is depicted in Fig. [Tl] We take care of the broken symmetry in the 
various steps of the derivation by setting the option symmetry-> broken in doRGE: 



Since we do not take into account the flow of the six-boson interaction, we do not use the option 
syimnetry->broken in doRGE. At this point this option should only be used if one is also interested in a deriva- 
tion of the flow equation of the six-boson coupling. This is due to the fact that the flow of the quartic interaction 
then depends directly on the flow of the six-boson coupling, see Ref. | 22| for details. 
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Figure 11; (Truncated) Flow equation for the two-point function in the case of broken 0{N) symmetry. The circles 
correspond to the non-vanishing expectation value of the field. 



twoPointB = doRGE[actionONSTriiiic2, {phi, phi}, symmetry -> broken]; 

The output is used as input for getAE: 

twoPointAlg = 

integrateDeltas [ 
Plus (SQ getAE [twoPointB, {{phi, il, 0, i}, {phi, ±2, 0, j»]]; 

The computation of the left-hand side is now slightly more involved than in the case with intact 
0{N) symmetry. We find 



To project onto the flow of the renormalized mass k we multiply both sides with S^^6^^ and replace 
the non- vanishing fields by po: 

twoPointAlgProj = 
integrateDeltas [delta [type , i, 1] delta [type, j, 1] twoPointAlg] /. 
phi[0, i_] :> phiO[i] /. phi0[_]~2 :> 2 rhoO 

— > -((2 lambda2b'2 rhoO dR[ql"2, k])/(R[ql~2, k] + ql"2 Z[ql"2])~3) + 
(2 lambda2b"2 rhoO dim [type] dR[ql'2, k])/(R[ql~2, k] + 
ql"2 Z[ql-2])-3 + (lambda2b dR[ql~2, k])/( 
2 (R[ql-2, k] + ql-2 Z[ql-2])-2) - (lambda2b dim[type] dR[ql"2, k])/( 
2 (R[ql-2, k] + ql-2 Z[ql-2])-2) + ( 

18 lambda2b"2 rhoO dR[ql"2, k])/(2 lambda2b rhoO + R[ql"2, k] + 

ql~2 Z[ql'2])"3 - (3 lambda2b dR[ql~2, k])/( 
2 (2 lambda2b rhoO + R[ql'2, k] + ql~2 Z[ql"2])~2) 

Let us now apply the rules dimLessRulesON to obtain a dimensionless fiow equation. This yields 

twoPointAlgProjy = 
Simplify [#, 
Element [d , 

Integers]] & /(§ (anglelnt[d] dql ql~(d - 1) twoPointAlgProj //. 
dimLessRulesON // Expand)// Simplify 




= ( - 6'^r]Zk{-n) + Zlk^-Ui{x)cPi{y){4^ -d-2r, + dt)X2 




(40) 



— > k-2 lambda2 y'(d/2) dtr [y] v[d] Z[k"2 y] "2 
(l/(y-2 (r[y] + Z[k-2 y])-2) 
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- dim [type] /(y-2 (r[y] + Z[k-2 y])-2) 

- 3/(y r[y] + (2 kappa lambda2 + y) Z[k"2 y])"2 

+ 4 kappa lambda2 Z[k~2 y] (-(l/(y'3 (r[y] + Z[k'2 y])~3)) 
+ dim [type] /(y'3 (r[y] + Z[k'2 y])~3) 
+ 9/(y r[y] + (2 kappa lambda2 + y) Z[k~2 y])~3)) 

On the left-hand side we project with S^^6^^: 



4>=<Po 



: (^-vZk{-n) + Zkk^X2{2^r]-d^dt)>i+ 

+ 2Zkk^K{4: -d-2r] + dt)\2)s{x - y). (41) 



In the point-Uke limit this reduces to 

dcj)\G)5(j)3 [0) 4>=<t'o 
= {Zkk^X2{2 -T]-d- dt)K + 2Zkk^K{A - d - 2?? + a^Aa) j d'^x. (42) 

Note that we have switched to momentum space to obtain this expression. The integral over d'^x 
is associated with the momentum conserving 5 function and cancels against the right-hand side, 
see also comments above. 

We observe that terms proportional to 2k^K appear on the left-hand side. These terms corre- 
spond to the flow equation of the quartic coupling. Inserting the RG flow equation of the quartic 
coupling into the expression above, we obtain the flow equation for the renormalized and dimen- 
sionless vacuum expectation value k: 

dtK^i2-7j^d)K~v, r dyy'^/'Zk ({l-N) ""'[^l 

Jo ^ V^{r{y) + ZkY 



3/(y) 



{yr{y) + {2k\2 +y)Zk) 



(43) 



Finally we have to specify a regulator function. To demonstrate how this can be done within 
DoFun, we consider the flow of the quartic coupling. In the following we choose an optimized 
regulator function Rk — (fc^ — q^)0{k^ — q^), see Refs. [TH l 11191 1120j . The momentum integrations 
can then be performed analytically with the aid of Mathematica. For this purpose, it is convenient 
to define the rules 

regRulesZ = {r[y] :> (1/y - 1) UnitStep[l - y] , 
dtr[y] :> 2/y UnitStep[l - y] , Z [_] :> 1}; 

Applying them to f ourPointT2Algy, i. e. to the integral appearing in the flow equation of the 
four-point function calculated above, we find: 

Integrate [ 

f ourPointT2Algy/(-3 k' (4 - d) ) /. regRulesZ, {y, 0, \ [Infinity] } , 
Assumptions -> {d > 0}] 



— > -((8 lambda2"2 (1 - 9/(1 + 2 kappa lambda2)~3 - dim[type]) 
v[d])/d) 

In standard notation the complete flow equation reads then 

= id-2 + 2,)X2 - 8A^^ ((1 - N) - . (44) 

Note that the same procedure can be applied to the expressions of all flow equations derived above. 
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5.2. Gross-Neveu-model 

In the previous section we have discussed a purely bosonic theory. In this section we turn to 
a purely fermionic formulation of the Gross-Neveu model. This model allows to study dynamical 
chiral symmetry breaking as driven by fermion fluctuations. In fact, the finite-temperature phase 
diagram of the Gross-Neveu model in d < 4 space-time dimensions has drawn a lot of attention 
in recent years |113l 1114] . Here, we only aim at a study of a purely fermionic formulation of this 
theory at vanishing temperature which is sufficient for a first non-trivial check of DoFun. 

The ansatz for the effective action we are going to employ is 

d d d 

+ ^^^^^('?i)^'('?2)^V;,(g3)V',(-9i-g3-g2). (45) 

We first define the fields psi and psib as fermion and anti- fermion, respectively, with a Dirac and 
a flavor index, 

defineFieldsSpecific [{{psi [momentum, dirac, flavor], 
psib [momentum, dirac, flavor]}}]; 

and then the action: 

actionGN = convert Action [ 
diracM[q, a, b] (Z[q"2] + Rs[q'2, k] ) op [psib [q, a, i] , 

psi[q, b, i]] + 
gb/ (2 Nf) op [psib [ql, a, i] , psi[q2, a, i] , psib[q3, b, j] , 

psi[q3 + ql - q2, b, j]]] ; 

In this expression Dirac 7 matrices are represented by the function diracM where the first argument 
denotes the momentum contracted with the 7-matrix and the other two represent the indices of 
the matrix. The wave- function renormalization is denoted by Z and the regulator function by Rs. 
The fiows of the two- and four-point functions can be obtained frorrp^ 

twoPoint = doRGE [actionGN, {psib, psi}]; 
fourPoint = doRGE [actionGN, 
{{psib, il}, {psib, 12}, {psi, 13}, {psi, 14}}]; 

The results are depicted in Fig. [T2| 

We now define the Feynman rules. From the action we obtain the algebraic expressions for the 
two- and the four-point functions 

FR2Polnt = getFR [actionGN, {psib[pl, a, 1], psi[p2, b, j]}] 

— > delta[f lavor , 1, j] deltam[-pl + p2] dlracM[pl, a, b] 
Rs[pl~2, k] + delta[f lavor, 1, j] deltam[-pl + p2] 
diracM[pl, a, b] Z[pl~2] 

and 

FR4Point = getFR [actionGN, {psib[pl, a, 1], psib[p2, b, j] , 
psi[p3, c, m] , psi[p4, d, 1]}] 

— > -((gb delta [dirac, a, d] delta [dirac, b, c] 
delta[f lavor , 1, m] delta [flavor, j, 1] 
deltam[-pl - p2 + p3 + p4])/Nf) + 
(gb delta [dirac, a, c] delta [dirac, b, d] 
delta[f lavor , i, 1] delta [flavor, j, m] 
deltam[-pl - p2 + p3 + p4])/Nf 



^^Note that for illustration purposes we directly use the physical instead of the symbolic action. 
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Figure 12: Flows of the two- and the four-point functions of the Gross-Neveu model. 



Wc use the expression for the four-point function directly and infer the form of the propagator 

from that of the two-point function: 

P[psi[pl_, a_, i_] , psib[p2_, b_, j_] , explicit -> True] := 
delta[f lavor , i, j] diracM[pl, a, b]/ 
pl~2/(Z[pl'2] + Rs[pl~2, k]); 
V[psib[pl_, a_, i_] , psib[p2_, b_, j_] , psi[p3_, c_, in_] , 
psi[p4_, d_, 1_] , explicit -> True] = 
FR4Point/deltaiii[-pl - p2 + p3 + p4] // Expand; 

Recall that DoFun implicitly assumes a momentum conserving i5 distribution to be present in its 
output. Furthermore we need to define the regulator insertion which we choose to be diagonal in 
flavor space: 

dR[psib[pl_, a_, i_] , psi[p2_, b_, j_] , explicit -> True] := 
diracM[pl, a, b] delta[f lavor , i, j] dRs[pl"2, k] ; 

Before we derive the algebraic form of the equations let us introduce dummy index names for 
Dirac (a, b, and so on) and flavor indices (i, j, and so on): 

addlndices [{{flavor , {i, j, k, 1, m, n}-}, 
{dirac, {a, b, c, d, e, f, g, h}}}] ; 

We also define some rules for calculations and convenient replacements later: 

dimRules = {dim [flavor] :> Nf, dim [dirac] :> dgamma}; 
diracRules = {diracM[-p_, a_, b_] :> -diracM[p, a, b] , 
diracM[p_, a_ , a_] :> 0, 

diracM[q_, a_, b_] diracM[q_, b_, c_] :> q~2 delta [dirac, a, c] , 
diracM[q_, b_, a_] diracM[q_, c_, b_] :> q"2 delta [dirac, a, c] , 
diracM[q_, a_, b_] "2 :> q"2 delta [dirac, a, a]}; 

The latter set of rules is useful for the evaluation of traces of products of Dirac matrices. 

Now we have all the ingredients at hand which are required to derive the algebraic expression 
for the two-point function with getAE: 

twoPointAlg = Plus 0® getAE [twoPoint, 
{{psi, i2, pi, b, j}, {psib, il, p2, a, i>}] // Expand; 

This is the flow equation of -rrb-, — xittt — 

^ (5V'°(pi)'5'0J(P2) 

For simplicity, we restrict ourselves to the point-like limit, i. e. we set the external momenta to 
zero. To simplify the expression twoPointAlg, we repeatedly apply the function integrateDeltas 
as well as the rules diracRules until the output does not change any further: 
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FixedPoint [integrateDeltas [#] /. diracRules &, 
twoPointAlg /. {pi -> 0, p2 -> 0}] 



— > (gb delta[f lavor , i, j] diracM[ql, a, b] dRs[ql'2, k])/ 
(Nf ql-2 (Rs[ql-2, k] + Z[ql-2])-2) 

Translated into the standard textbook notation, this expression reads 



As expected, this expression vanishes in the point-hke Umit since the integrand is anti-symmetric 
in the four- momentum q^. 

The flow equation of the four-point function can be derived along the lines of the derivation of 
the two-point function. For convenience, however, we set the external momenta to zero right from 
the beginning: 

fourPointAlg = Plus 0® getAE [f ourPoint , 

{{psib, 11, 0, a, 1}, {pslb, 12, 0, b, j}, 
{psl, 13, 0, c, 1}, {psl, i4, 0, d, m»] ; 

Applying the same simplification procedure as above yields 

f ourPointAlgSlm=Simplif y [FixedPoint [integrateDeltas [#] / . 
diracRules &, fourPointAlg]] 

— > (l/(Nf-2 ql-2 (Rs[ql-2, k] + Z [ql-2] ) "3) )2 gb-2 

(-2 + dgamma Nf) (delta[dlrac, a, d] delta[dirac, b, c] 
delta[f lavor , i, m] delta [flavor, j , 1] - 
delta[dirac, a, c] delta[dirac, b, d] delta[f lavor , i, 1] 
delta[f lavor , j, m] ) dRs[ql~2, k] 

This is the right-hand side of the flow equation of the four-point function g^ag^b/^a^d ■ There are 

only two terms left. Thus, contributions from the three different diagrams must have canceled each 
other. Since we would like to compute the flow of the coupling g, we have to project the output 
f ourPointAlgSim onto our ansatz Eq. (45 I. To this end we contract the output f ourPoint AlgSim 
with the tensor-structure of the four-fermion interaction given by 

V[psib[pl, a, i] , psib[p2, b, j] , 

psi[p3, c, 1], psi[p4, d, m] , explicit -> True]/gb 

This yields the flow equation for g: 

f ourPointAlgSimProj = FixedPoint [integrateDeltas [#] /. 
diracRules &, V[psib[pl, a, i] , psib[p2, b, j] , 

psi[p3, c, 1], psi[p4, d, m] , explicit -> True]/gb 
f ourPointAlgSim] /. dimRules; 
Ihs = integrateDeltas [ 

V[psib[pl, a, i] , psib[p2, b, j] , psi[p3, c, 1], 

psi[p4, d, m] , explicit -> True] 
V[psib[pl, a, i] , psib[p2, b, j] , psi[p3, c, 1], 
psi[p4, d, m] , explicit -> True]/gb~2] /. dimRules; 

Ihs corresponds to the projection of f ourPointAlgSim on the left-hand side. Thus, the flow 
equation for the four-fermion interaction g reduces to the following simple expression 

beta = f OurPointAlgSimProj /Ihs // Simplify 

— > -((2 gb"2 (-2 + dgamma Nf) dRs[ql~2, k])/ 
(Nf ql-2 (Rs[ql-2, k] + Z[ql-2])-3)) 

Let us rewrite this expression with the aid of the following replacement rules: 
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dimLessRulesGN = {Rs[ql~2, k] :> r[y] Z[ql~2], gb -> g k~ (2 - d) , 
dRs[ql~2, k] : > -2 y Z[ql~2] r' [y] , ql :> Sqrt [y] k, 
dql :> k-2/2/ql}; 

Here, r [y] denotes the dimensionless regulator function, g the dimensionless couphng constant 
and y=ql~2/k~2. dql is the integral measure. Theses rules allow us to obtain a flow equation in 
terms of dimensionless quantities. The integral over the angular part is 

anglelnt[d_] :=4v[d]; 

where v[d] represents Vd — {n'^^'^2'^~^^r{d/2))~^ . Next, we rewrite the flow equation in terms of 
dimensionless quantities: 

betaDimLess = 

angle Int [d] dql ql"(d - 1) beta //. dimLessRulesGN// Expand; 
Simplify [betaDimLess , Assumptions -> Element [d, Integers]] 

— > (8 g~2 k'(2 - d) (-2 + dgamma Nf) y~(-l + d/2) v[d] r'(y))/ 
(Nf (1 + r[y])-3 Z[k-2 y]-2) 

Translated into standard textbook notation this expression reads 

{d^Nf-2)k — - / dy—-j-^—--^. (47) 



Nf Jo "(l + r(y))3Z(fc2y)2 

To obtain the final expression for the flow equation we also have to manipulate the left-hand 
side. It reads 

' -An- ' 



{{2 - d - 2j]^)g + dtg) k 



2-d 



X 



X ^ {5u5,^5'^^5''' - 5,^5,i5'^H'-) {2n)H{p + q + r + s) (48) 
iV/ 

where 77^, = ~dt \nZ is the anomalous dimension associated with the field i/j. Note the minus sign 
in front of the derivatives which appears due to our conventions, see Eq. ([s]). In the point-like limit 
the projection yields 

((2 - d - 277v,)<? + dtg) A:2-'^(27r)'^5(0). (49) 
Combining the left- and right-hand sides we obtain the flow equation for g: 

d,g^id-2 + 2,,)g + (d, - - j 8.. dy (^^^(^5^. (50) 

In the large iV"/-limit this reduces to the result found in Ref. [121| . 

Of course, the six-point function can be also derived with DoFun. From the symbolic expression 

sixPoint = doRGE [actionCNSymbolic , 
{psib, psib, psib, psi, psi, psi}] ; 

(which yields 45 terms as can be checked with countTermsQsixPoint) we obtain the algebraic 
expression: 

sixPointAlg = Plus 00 getAE [sixPoint , 

{{psib, il, 0, al, jl}, {psib, i2, 0, a2, j2}, {psib, 13, 0, a3, j3>, 
{psi, 14, 0, a4, j4}, {psi, 15, 0, a5, j5}, {psi, 16, 0, a6, j6»] ; 

Simplifying the result reveals that each term contains an odd number of ^ and thus the flow of the 
six-point function vanishes upon integration in the point-like limit as expected. 

To illustrate the flexibility of the developed formalism we would like to mention how easy it is to 
extend our study to non-vanishing external momenta. With DoFun, we only need to change a few 
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arguments and define some additional rules for the handling of the Dirac algebra. We demonstrate 
this for the flow equation of the four-point function. In short, we perform the following steps: 
Derive the algebraic RG equation, 

f ourPointAlgFinMom = 
Plus m getAE[ 

fourPoint, {-[psib, il, pi, a, i}, {psib, 12, p2, b, j}, 
{psl, 13, p3, c, 1}, {psl, 14, p4, d, m}}] ; 

project it using 

f ourPolntAlgProj = 
FlxedPolnt [IntegrateDeltas [#] /. dlracRules &, 

V[pslb[pl, a, 1], pslb[p2, b, j] , psl[p3, c, 1], psl[p4, d, m] , 
explicit -> True]/gb f ourPolntAlgFlnMom] /. dlmRules; 
betaPlnMom = f ourPolntAlgFlnMom/lhs ; 

and then make it dimensionless, 

betaPlnMomDlmLess = 
anglelnt[d] dql ql~(d - 1) betaPlnMom //. dlmLessRulesGN // Expand; 

For simplicity, we finally present the result for the large- A'f limit: 

betaDlmLessLargeNf = Limit [betaDlmLess , Nf -> \ [Infinity] ] ; 
Simplify [betaDlmLessLargeNf , Assumptions -> Element [d. Integers]] 

--> -(4 g-2 k-(2 - d) y-(-l + d/2) 

dlracM[k Sqrt [y] , C$180535, f $180535] v[ 
d] ((pi + p3 + k Sqrt [y]) "2 dlracM[-p2 - p3 - k Sqrt [y] , 
f$180535, 

C$180535] (Rs[(pl + p3 + k Sqrt[y])~2, k] + 
Z[(pl + p3 + k Sqrt[y])~2]) + (p2 + p3 + 
k Sqrt[y])~2 dlracM[-pl - p3 - k Sqrt [y] , f $180535, 
C$180535] (Rs[(p2 + p3 + k Sqrt[y])~2, k] + 
Z[(p2 + p3 + k Sqrt[y])~2])) Derivative [1] [r] [ 
y])/((pl + p3 + k Sqrt[y])-2 (p2 + p3 + k Sqrt[y])-2 (1 + 
r[y])~2 (Rs[(pl + p3 + k Sqrt[y])~2, k] + 

Z[(pl + p3 + k Sqrt[y])~2]) (Rs[(p2 + p3 + k Sqrt [y]) "2, k] + 
Z[(p2 + p3 + k Sqrt [y]) -2]) Z[k-2 y] ) 

This result can be further simplified and used in numerical calculations. 
6. Summary 

In this paper we have presented DoFun, a Mathematica application which allows to derive both 
Dyson-Schwinger equations and functional RG equations starting from a given action. DoFun is 
based on DoDSE, see Ref. [98], which was limited to the derivation of DSEs in symbolic form. Do- 
Fun goes beyond the symbolic form and provides explicit expressions for the integrals. Apart from 
the functions used for the actual derivation of the equations, we have included several additional 
helpful tools, e. g. for dealing with the Kronecker delta and for the derivation of Feynman rules. 

We have demonstrated the usage of DoFun by means of a scalar 0{N) field theory and the 
Gross-Neveu model in a purely fermionic description. In particular, 0{N) symmetric scalar field 
theories have been studied in great detail in the literature, see e. g. Refs. [22H26] . Although 
these theories are well suited to demonstrate the usage of DoFun, our main goal was to provide a 
tool which facilitates the derivation of functional equations for even more involved theories, such 
as QED or QCD, where the tensor structure leads to additional complications. Furthermore, the 
numbers of terms grow considerably in gauge theories when higher n-point functions are taken into 
account. Therefore we hope that DoFun proves to be a valuable tool in particular for functional 
studies of gauge theories and helps to push them to new limits. 

We would like to encourage users of DoFun to actively communicate their experiences with it in 
order to help us to further improve this application. In particular, bug reports are most welcome. 
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Appendix A. How to include Grassmann fields 

Implementing the anti-commuting nature of fermionic fields into a Mathematica program like 
DoFun requires the use of a non-commuting multiplication operator. This was one reason for 



introducing the op-function as described in section 3.2 It allows to properly take into account 
the effects of changing the order of the fields. Internally all anti-commuting fields are called either 
fermion or antiFermion. This also includes, e. g., ghost fields, which are not fermions, but 
anti-commute as well. The corresponding test functions are f ermionQ and antiFermionQ, while 
grassmgoinQ yields true in either case. 

DoFun recognizes fermions via the definition of the two-point functions: If any two-point 
function in the action has two fields with different names, DoFun considers the first one as a fermion 
field and the second one as its anti- field. In case of complex fields or a mixing of bosonic fields at the 
level of two-point functions this can be suppressed with the option specif icFieldDefinitions, 
see section |3.1| DoFun also allows to define all fields as bosons or fermions without specifying 
two-point functions. An example for a fermionic action is given by the Gross-Neveu modep] 

actionGNSymbolic={{psi , psib}, {psib, psib, psi, psi}}; 

Here, the canonical order was used where anti- fermions (psib for ip) appear always to the left of 
fermions (psi) in n-point functions. This is due to the fact that we employ the following definitions 
for Grassmann derivatives: 

(5 5 S ~S , 

— M := A/ — , -^M := -^M. (A.l) 

dip dip oip dip 

Consequently, derivatives with respect to fermions are assumed to act from the right whereas 
derivatives with respect to anti-fermions are assumed to act from the left. For fermions, the 
differentiation rules of Eq. ([t]) have to be adapted as follows: 

^_ T^jk _ ^_ ( ^ ^ \ _ nJ™ j^imn j~,nk 



To alleviate the readability of expressions, propagators are defined with fermions and anti- 
fermions exchanged, i. e. a propagator Df^{p'^) is in DoFun correctly represented by P[{psi,i}, 
{psib.j}]. As a consequence, the connecting fields of a vertex and a propagator are the same 
even for fermions, whereas normally one would have a fermion and its anti-fermion. For example, 
in the two-point DSE of the Gross-Neveu model we encounter the following tadpole-like diagram 
where the leg {psib.rl} of the vertex is attached to the leg {psib.rl} of the propagator: 

op[S[{psib, i2}, {psib, rl}, {psi, si}, {psi, il>] , 
P[{psi, si}, {psib, rl}]] 



'In this section i/) and i/) correspond to Grassmann fields, while <I> is a bosonic field. 
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We choose this notation to make it easier to connect the propagators and vertices visually when 
seeing this output. 

Using left- and right-derivatives also has consequences for the definition of vertices. For instance, 
a four-fermion interaction has the form 



-i.2 <_2 



Wphys •,'Pphys 



Thus, the indices in fVt', , do not reflect the order in which the derivatives are performed, but 
rather the order in which the derivatives appear. In this example this means that the differentiation 
with respect to {ipk) has to be performed before that with respect to ipi {tpi). This convention 
is employed for bare and dressed vertices, S and V, respectively, and also in getFR. To derive the 
expression for the bare four-point function 5*^^^^ from a given action one uses 

getFR[action, {psib[pl,i], psib[p2,j], psi[p3,k], psi[p4,l]}] 

The required reordering of anti-fermions is done automatically by getFR. 

An exception to that rule for fermion ordering are the derivative arguments in doDSE. Here, 
the order corresponds exactly to the order in which the derivatives are performed. This is due to 
the fact that the form of a DSE depends on which field is attached to the bare vertex. For the 
ghost-gluon vertex of Landau gauge Yang-Mills theory this is explicitly demonstrated in ref. [98] . 
For example, the DSE of the four-point function F^^-^^^ can be obtained from 

doDSE [action, {{psib, j}, {psib, i}, {psi, k}, {psi, 1}}] 

but also from 

doDSE [action, {{psi, k}, {psi, 1}, {psib, j}, {psib, i}}] 

Other variations also exist. The important point is that the first derivative has to be taken either 
with respect to "(Afc or ■0^. For RGBs the same rule applies. 

Now we turn to the minus signs arising from the anti-commutativity of fermions. As we use 
superfields which may have bosonic as well as fermionic entries, some care is required in the 
derivation. For DSEs the minus signs for fermionic loops arise immediately from the ordering of 
the fields at the end of the calculation. An example in ref. [55] nicely demonstrates this. However, 
the currently employed algorithm is not infallible: There are known examples with a wrong sign, 
but the lowest diagram affected is a two-loop diagram of a three-point function. In this case the 
sign needs to be corrected by hand. 

For RGEs we employ a completely different algorithm and up to now we are not aware of any 
shortcomings. In contrast to DSEs the order for the field derivatives does not affect the result. 
Therefore we are free to put them into canonical order: anti-fermions always appear to the left of 
fermions to the left of bosons. The order between fields of the same 'type' is not changed. Hence 
doRGE always uses this canonical order to derive RGEs. This ordering is in agreement with our 



choice given in Eq. ( 13 1. For the differentiation itself we ignore any minus signs that would arise by 
dragging one Grassmann derivative to the left or to the right of another one. The minus signs from 
passing two Grassmann fields of the external derivatives are taken into account at the end. As far 
as the internal fields are concerned, the nature of the fields (bosonic or Grassmann- valued) is of no 
importance since after setting the sources to zero all propagators, vertices and regulator insertions 
have zero Grassmann number and can be freely commuted. In fact, we find that all minus signs 
that would arise if we took into account the Grassmann numbers of internal fields cancel each 
other. Thus, the only signs we have to worry about stem from interchanging the external field 
derivatives. Those are corrected in a last step by comparing their order with that of the original 
canonically ordered derivatives. 



Appendix B. Version overview of DoFun 

The predecessor of DoFun is DoDSE, see Ref. |98i. Subsequently smaller updates have been 
made available. With version 2.0 the package DoDSE became part of the application DoFun 
together with the new packages DoAE and DoFR. The current version is 2.0. 

The following list gives a short overview over all publicly available versions: 
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• DoDSE 1.0 (Aug. 15, 2008): first publicly available version 

• updates DoDSE 1.1, 1.2, 1.2.1, 1.2.2, 1.2.3 (Feb. 9, 2010) 

• DoFun 2.0 (Feb. 25, 2011): contains DoDSERGE (formerly DoDSE), DoAE and DoFR 
Further updates will be made available at 

http: //theorie . ikp.physik. tu-darmstadt . de/~mqh/DoFun[ 



If an Internet connection is available, DoFun will automatically notify the user about available 
updates. 
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